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Starting from a simple classical framework and employing some stochastic concepts, the basic 
ingredients of the quantum formalism are recovered. It has been shown that the traditional ax- 
iomatic structure of quantum mechanics can be rebuilt, so that the quantum mechanical framework 
resembles to a large extent that of the classical statistical mechanics and hydrodynamics. The 
main assumption used here is the existence of a random irrotational component in the classical mo- 
mentum. Various basic elements of the quantum formalism (calculation of expectation values, the 
Heisenberg uncertainty principle, the correspondence principle) are recovered by applying traditional 
techniques, borrowed from classical statistical mechanics. 
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I. INTRODUCTION 

One of the most intriguing achievement of the 20-th 
century physics is the foundation of quantum mechanics 
and its basic tool the Schrodinger equation. It is one 
of the most studied equation in contemporary physics 
both from mathematical point of view, as well as from 
the perspective of the enormous number of its important 
applications. Various links have been proposed between 
classical and quantum picture, in order to overcome the 
main difficulties due to the difference of formalisms, and 
to better understand their interplay and possible connec- 
tion. 

A formal approach is represented by the Wigner-Weyl 
transformation , based on the reformulation of quan- 
tum mechanics into phase space by means of a continu- 
ous map. Although having some drawbacks, such as the 
fact that the Wigner distribution function is not strictly 
positive, this approach has proved an important tool in 
many areas of quantum physics and chemistry. In partic- 
ular, the semi-classical limit of quantum mechanics can 
be recovered by using the Wigner-Weyl transformation. 

While the Wigner-Weyl approach can be considered 
as an extension of quantum mechanics in a classical 
domain (phase space), the Koopman-von Neumann ap- 
proach H H Q is a reformulation of classical mechanics 
into the Hilbert-space language. This quantum-like the- 
ory, although being always classical has the advantage 
that it can be directly compared to quantum mechanics, 
at least in terms of the underlying formalism. 

The two guidelines of development mentioned above 
are actually the two sides of the same subject, namely 
to find a formal mathematical language capable to han- 
dle both classical and quantum mechanics. There ex- 
ists a third attempt fundamentally different compared to 
those, which is known as the Bohm-Fenyes-Nelson ap- 
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This approach is centered 
on the physical footings and interpretation of quantum 
mechanics. The basic idea is that quantum phenomena 
are conditioned by stochastic effects, which take place 
in a classical framework, so that the notion of equations 
of motion and hence the notion of trajectory (although 
a random one) remains valid. The Planck's constant h 
plays now the role of a measure of the strength of the 
stochastic effects. This is an interesting point of view 
exhibiting a number of remarkable properties and links 
to the theory of deformation quantization. Many arti- 
cles H2lhJ,L4, 151 are devoted to the study of the rela- 
tion between the phase space representation of quantum 
states (quasi-distributions) in the deformation quantiza- 
tion method, and the Bohm distributions in phase space 
in the framework of the Bohm-Fenyes-Nelson approach. 

Recently, it has been shown ^(| that an exact uncer- 
tainty principle can be formulated, which provides the 
key argument in the transition from the dynamical de- 
scription of a classical ensemble to that of a quantum 
ensemble. Another interesting derivation of the equa- 
tions of nonrelativistic quantum mechanics is based on 
the use of the principle of minimum Fisher information 

m 

The point of view pursued in the present paper has 
been formulated previously in a different context |l8| . Its 
cornerstone is the assumption that the particle velocity, 
being an infinitesimal quantity (derivative of the position 
with respect to time), is random and consists of a mean 
part and an irrotational fluctuation. The mean part is the 
actual classical velocity (momentum), while the fluctua- 
tion represents a measure of the uncertainty with which 
it can be determined. Since the fluctuation is irrota- 
tional, it does not affect the averaged Hamilton equations 
of motion, however it yields an additional second-order 
correlation term in a picture, where an ideal particle lo- 
calization is not valid. Moreover, as it will become clear 
in the sequel, the strength of the momentum fluctuation 
depends functionally on the measure of particle dereal- 
ization, namely on the density distribution g(x;t). 
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The paper is organized as follows. In the next Section 
a brief parallel between the Hamiltonian and Liouvillean 
description is outlined. In Section III, a mapping of the 
classical kinetic balance equations onto a Schrodinger 
equation is derived. The second-order correlation ten- 
sor of the random, irrotational velocity field is derived 
in Section IV. Further, it is shown that the hydrody- 
namic equations thus obtained, are in fact the Madelung 
equations 0, known to be formally equivalent to the 
Schrodinger equation. Section V presents a further elab- 
oration aimed to recover some basic rules of the quantum 
formalism. Finally, Section VI is devoted to discussion 
and conclusions. 



In addition, the function <S(x; t) satisfies the Hamilton- 
Jacobi equation and represents a family of classical tra- 
jectories. Evidently, the density distribution g(x; t) is a 
new element in general. It can be regarded as a gen- 
eralization of the delta-distribution in the Hamiltonian 
description, and subsequently as a measure of particle 
derealization. 

Let us define the characteristic function 

£(x, s; t) = J dp/(x, p; t) exp fjS ■ p^j , (II.6) 
with a Fourier inverse 



II. LIOUVILLEAN DESCRIPTION 

To start with, we consider a n-dimensional Hamilto- 
nian system, whose dynamics is governed by a smooth 
Hamiltonian function H(x.,p;t). Here x and p are n- 
dimensional vectors. It is well-known that the descrip- 
tion of the dynamical system in terms of equations of 
motion (Hamilton's equations) is formally equivalent to 
the description in terms of a phase space density specified 
by a distribution function /(x, p;t). The latter satisfies 
the Liouville equation 



where 



{F,G} 



>p / dF_ dG_ _ dG dF \ 
£-J V dx k dp k dx k dp k J 



(III) 



(VsF) • (V P G) - (V a G) • (V P F), 



(IL2) 



is the standard Poisson bracket. This equivalence is how- 
ever subtle. Actually, it is a formal mathematical prop- 
erty, which is evidently incomplete in a physical sense. 
Note for example that, unlike the Hamilton's equations 
of motion 



(II.3) 



where x = x(t) and p = p(£), in equation ljll.ll) x and 
p are independent variables and they do not depend on 
t. Moreover, the Hamilton equations correspond to the 
ideal case of a perfectly localized particle 

/(x, p; f) = 6[x - X.(t)]8\p - P(t)], (II.4) 

where (X(i),P(i)) represents their explicit solution. Ex- 
pression ljll.41) also satisfies the Liouville equation, how- 
ever in addition, it admits a solution of the form (|IV.15(I . 
The latter implies that for a given value of x, the clas- 
sical momentum p is uniquely determined according to 
the expression [see equation i|IV.34[l ] 



/(x,p;t) 



dsQ(x, s; t) exp 



A 



s • p 



(II.7) 

where the variable s is chosen such that to have the same 
dimension as the coordinate x, and A is a formal param- 
eter with dimension of action. It is introduced in order 
to make the argument under the exponent dimensionless. 
In fact, equation (|II.6(I [or equivalently (|II.7ll ] represents 
the definition of the Wigner-Weyl transformation • For 
the time being, we avoid any reference to the latter, apart 
from formal similarity and coincidence. In Section V, we 
will establish a firm link between the characteristic func- 
tional and the Wigner-Weyl representation. 

Before we proceed further, let us point out a couple 
of important features of the characteristic function. It is 
known that in the limit s — > 0, the characteristic function 
yields the moments of the distribution function /(x, p; t) 
after a proper marginalization with respect to the conju- 
gate momentum variable p is performed. For example, 



hm£(x,s;*) = 



where m is the particle mass, and 



(II. 



Kf V N 

N= hm — , 

N,V^oo V 



(II.9) 



p = V x <S(x;t). 



(II.5) 



is the particle number density in the thermodynamic 
limit. Here N is the total number of particles in the 
system and V is the volume occupied by the system. Fur- 
ther, we have 



- lim V s £(x, s; t) = -U(x; t)v(x; t), (11.10) 

l s^O J\l 



and so on. 

It is therefore very instructive to derive an equation 
for the characteristic function and try to manipulate it 
in a suitable manner. 
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III. ANALYSIS OF THE LIOUVILLE 
EQUATION AND MAPPING ONTO A 
SCHRODINGER EQUATION 

In this Section, we will work out in detail the simplest 
case, where 



£r(x,p;t) = f- + tf(x;t). (III.l) 
2m 



For each term in the Liouville equation ijll.lfl . we have 
subsequently 



df i r.dg ( i 



(III.2) 



V,/ = J dsV s • V X G exp (-is • p) , 

(III.3) 



Therefore, the characteristic function can be written ac- 
cording to the relation 



C/(x, s; t) = F(x, s; t) exp 



-G(x,s;i) 



where 



(III. 12) 



F(x, s;t) = Ui (x + ~;t)Kz(x ~ |; *), (HI.13) 
G(x, s; t) = S(x + |; i) - s(x - |; t). (IIL14) 



From equation (|111.14|l it becomes clear that the phase 
G(x, s;t) is an odd function of the variable s. Let us 
now substitute the characteristic function represented by 
(|III.12(1 into equation (|III.5(1 and separate the real and 
imaginary part. We obtain 



ds(s ■ V X U)G exp 



-s • p 



A(27rA)" J y ' 1 V A 

(III.4) 

Combining all the terms, the sought-for equation can be 
written the in form 



at m A 



(III.5) 

We would like now to diagonalize the differential opera- 
tor, encountered in the second term of equation (|III.5(I . 
This is achieved by a simple linear change of variables 



x 1= x+-, 



x 2 = x . 

2 2 



Taking into account the identities 
V 2; =Vi+V 2 , V s =-(Vi 



(III.6) 



(III.7) 



where V1.2 denotes the differential operator taken with 
respect to the variables Xi j2 respectively, we note that 
equation l|III.5|) can be rewritten as 

% - ^(vf - vl)e + i(s • vjj)q = 0. (111.8) 



The last equation suggests the ansatz 

£(xi,x 2 ;t) = * 1 (x 1 ;t)* 2 (x 2 ;i). 



(III.9) 



In addition, the yet unknown complex valued functions 
(fc = 1, 2) can be represented in the form 



* fc (x fe ;i) = 7£ fe (x fc ;t)exp 



A 



1,2. 



(111.10) 

Since <?(x, s = 0;t) must be real, it follows immediately 
that 



S(x;i) =<S!(x;t) = -S 2 (x;t). 



(III.ll) 



+ V x ■ (FW S G) + (V S F) ■ (W X G) = 0, (111.15) 



dC F A 2 
F— + ~(V.G) • (V K G) = -Fa ■ V X U + — V s • V^F. 
at m m 

(111.16) 

In what follows, we will analyze equations (|III.15|I and 
(|III.16(1 in zero and first order with respect to the s- 
variable. The zero-order reads as 



dt 



-v x -(f^ 

m V 



V.5 =0, 



K 2 V z fti - ^iV x ^ 2 = 0, 



(111.17) 



(111.18) 



where F(°'(x;f) = fti(x;i)ft 2 (x;i). Equation (|III.18|) 
simply implies that 1Z\ and 7?. 2 must be equal up to a 
multiplicative constant that can be normalized. Thus, 
without loss of generality, we can write 



K(x;t) =7£i(x;t) = 7£ 2 (x;t). 



(111.19) 



Therefore, the amplitude F(x, s; t) is an even function of 
the s-variable 

F(x, s; t) = k(x + l;t)n(x - |; t) , (111.20) 
In first order, we obtain 



VJ — + ; = -V X U H V- — ^— , . 

x \dt) x 2m x 2m x \ K J 

(111.21) 

which integrated once yields 

95 ■ (V ^ )2 - -tZ+^S?. (111.22) 



at 2™ 



2m K 
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The final step is to introduce the complex wave function 
tp(x; t) according to the de Broglie ansatz 



t) = 7?.(x; t) exp 



-S(x;t) 



(111.23) 



Combination of equations l|III.17|) and ljIII.22|) yields the 
following result 



\ 2 

2m 



Uip. 



(111.24) 



By identification of A with the Planck's constant H, equa- 
tion (|111.24|) transforms into the Schrodinger equation. 
In addition, equations (|III.17|I and l|III.22|) coincide with 
the system of equations, describing the properties of the 
Madelung fluid |19|. 

The above considerations can be repeated for the case 
of nonrelativistic motion of a spinless particle in electro- 
magnetic field, governed by the Hamiltonian 

ff(x, p; *) = J_[p - e A(x; t)f + e£/(x; t). (111.25) 
2m 

It is worthwhile to mention that the approach based on 
the characteristic function yields the Schrodinger equa- 
tion up to first order in the variable s. In this sense it can 
be considered as an infinitesimal mapping of the classical 
kinetic balance equations onto the Schrodinger equation 
as pointed out previously pol |2]| . However, a drawback 
of this method is evident since ansatz (|III.9|I restricts a 
possible class of states, while some classical distributions 
are not described with it. 



IV. IRROTATIONAL MOMENTUM 
FLUCTUATIONS 

Our basic assumption concerns the equitability of po- 
sition and momentum, which are obviously not on the 
same footing. We assume position to be a fundamen- 
tal variable, while momentum being proportional to the 
infinitesimal variation of position respective to an in- 
finitesimal variation of time cannot be determined ex- 
actly. The physical argument for such assumption is the 
following. If an object is perfectly localized, there is no 
reason for the impossibility to determine its velocity ac- 
curately. If however, a probability assignment in config- 
uration space strongly violating particle localization is at 
hand, there must be some uncertainty in the specification 
of the infinitesimal variation of the particle " position" in 
the course of time. 

Following instead of (|III.1|I . we consider a dynam- 
ical system described by the Hamiltonian 



Defining the new variable 

P = P + A, (IV.3) 
we can write the Hamilton equations of motion as follows 

dx P dP dA „ TT P „ t /TTr 
-77 = -, — = ---V x U--xV x xA.(WA) 
at m at at m 

Suppose now that A is irrotational, that is 

V x x A = 0, (IV.5) 

which also implies 

A = -V B $. (IV.6) 

With these observations, it follows that the Hamilton 
equations l|IV.4|) can be obtained from a new Hamilto- 
nian 

W (x,P;t) = ^+tf(x;t) + ^^, (iv.7) 

Next, we perform a polynomial marginalization if the 
distribution function f(x,P;t). This is done by initially 
multiplying the stochastic Liouville equation 



Of 



dt m 

lfcl r,k 



V,/ + F-V P / = 



(IV.8) 



by various powers P* 1 P£ 2 ■ ■ ■ P^J 1 an d then formally inte- 
grating over the momentum variable. Here, the random 
force F is given by the expression 



F = ^-V X U. 
dt 



(IV.9) 



The equations for the first two moments can be written 
as 



|f + V, • (gV) = 0, 



(IV.10) 



d(gV n ) | dU kn _ gFn 
m 



dt 



(IV. 11) 



where 



p(x; t) = mjV J dP/(x, P; t), (IV.12) 



A(x;f)] 2 + [/(x;t), (IV.l) p ( X ; i)V(x; t) = nJ dPP/(x, P; t), (IV.13) 



where A(x;t) is yet unspecified fluctuating part of the 
classical momentum with vanishing mean value 



(A(x;t)> =0. 



(IV.2) 



n fci (x;i) = — / dPP fe F ; /(x,P;i). (IV. 14) 

m J 
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It is a simple matter to verify that the Liouville equation 
(|1V.8|I possesses an exact solution of the form 

/(x,Pif) = ^5[P-mV(x;t)], (IV.15) 

which is also known as the classical Bohm distribution, 
usually interpreted as the classical limit of quantum pure 
states Substitution of the classical Bohm distribu- 
tion into equation HIV. 14(1 yields the expression 



Il kl =gV k V h 



(IV. 16) 



for the stress tensor Hjy. Hence the system <|1V.10[) 
and (|iV.ll|) represents an exact closure of hydrodynamic 
equations, fully equivalent to the Liouville equation. 

We now take into account the fact that the momentum 
P and hence the current velocity V consists of a classical 
mean part v (corresponding to p) and a fluctuation V 
(corresponding to A). Averaging equations (|IV.10|) and 
(|IV.llfl . we obtain 



8g 

at 



+ V x • (gv) = 0, 



dv n 

dt 



Vk 



dv n 

dx k 



1 dU 

m dx n 



gox k 



(IV. 17) 



(IV.18) 



where 



(VkVn 



(IV.19) 



Recall that it was initially assumed that the momen- 
tum fluctuations are irrotational. This also implies irro- 
tationality of the second term on the right-hand-side of 
equation l|iV.18(l . It can be written in the form 



8Z 

dx n 



0C h 



kn. 



dxi 



OR 

dx k ' 



(IV. 20) 



where 



R = lng. (IV.21) 
Taking curl of both sides of equation l|IV.20(l , we obtain 



eh, 



dx 



d fOC kn 8R\ 

-R hCfc„- = U, 

Ox k ox k J 



(IV.22) 



where as usual, e; mn denotes the fully antisymmetric 
third-rank unit tensor. Multiplication by ei qp and sum- 
mation on I in the last identity, yields a second order 
linear partial differential equation for the unknown cor- 
relation tensor Ckn 



d | dR\fdC kn dC km 

dx k dx k ) \ dx m dx n 



Ckn 



d 2 R 

dxkdx r , 



c d2R -n 

o ^ o. 

OX k OX n 



(IV.23) 



r 



Note that the left-hand-side of equation (|IV.23|) is an- 
tisymmetric with respect to the indices m and n. This 
restricts considerably the number of its solutions. On the 
other hand, equation ljiV.23() is a linear equation with re- 
spect to Ck n , so that its general solution can be written 
as a linear combination of particular solutions. Since the 
correlation tensor is symmetric, it can be represented in 
diagonal form. To find a particular solution, suppose 
that equation (jIV.23|l is written in a reference frame in 
which the correlation tensor is diagonal. Then, all el- 
ements must be equal, which is the only possibility for 
this particular solution [represented by the second term 
in equation (|IV.24I) ]. Another solution can be represented 
in the form of a Hessian matrix of a generic function. A 
simple and straightforward verification shows that these 
two particular solutions exhaust all possibilities, and Ckn 
can be written as follows 



Ci 



d 2 v 

dxkdx 7 



+ flSknF, 



(IV.24) 



where a and j3 are constant coefficients, and the scalar 
functions ^(x; t) and T(x; t) are some functions of x and 
t that remain to be determined. Direct substitution of 



the solution l|IV.24|) into equation ljiV.23|) shows that J- 
is an arbitrary function of R 



^(x;t) = F(R), 



(IV.25) 



while r is equal to R up to a multiplicative constant, 
which (without loss of generality) can be set equal to 
unity 

r(x;i) =R(x;t). (IV.26) 
Hence, equation I|1V.20|I can be rewritten as follows 



V X Z = aS7 x 



2aV x 



V 2 x R+±(V x Rf 



\7 X T + T\J X R 



V X T + TV X R. 



(IV.27) 



We would like now to show that the arbitrary con- 
stant a must be negative {a < 0). In analogy to the 
definitions (|iV.12|l - ijiV.14() . we can introduce the kinetic 
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energy density according to the relation 

£(x; t) = ^- f dPP 2 /(x, P; t) = (IV.28) 
Averaging equation (|IV.28|I . we obtain 

(£?(x;t)) = ^ + |TW, (IV.29) 

where TrC denotes the trace of the correlation tensor 
(|IV.19(I . The second term on the right-hand-side of equa- 
tion (|IV.29(1 represents the density of the internal energy, 
which is due to the fluctuating part of the current veloc- 
ity. The total internal energy is given by the expression 



1 



dxpTrC 



dx 



Q 

(IV.30) 

where integration by parts and taking into account of 
vanishing integrals has been performed in the second 
term on the right-hand-side. Since the total internal 
energy must be positive for any choice of the arbitrary 
function T (including T = 0), the free parameter a must 
be negative. Remarkably enough, the second term on the 
right-hand-side of equation (|iV.30|) is proportional to the 
Fisher information 22] 



X = / d 



(V^) 2 



(IV.31) 



The above equation (|IV.31|) represents a direct link be- 
tween quantum mechanics and Fisher information theory. 
It is remarkable that in the simplest case, where T = 0, 
the Fisher information measure X has a very transparent 
meaning as being simply a measure of the internal energy 
of the quantum system. 

First of all, we would like to explore the simplest case, 
where T — 0. Identifying the parameter a as 



a 



Am 2 ' 



equation l|IV.18|) can be written accordingly 

— + v • V x v = + — V, — j=- 

ot m 2m z \ y/g 



1 



If we further define 

v = -V X S, 
m 

the de Broglie ansatz (|III.23() with 

A = h, K= y/g, 

yields immediately the Schrodinger equation 



ih 



9-0 
df 



2m 



(IV.32) 



(IV.33) 



(IV.34) 



(IV.35) 



(IV.36) 



Another simple but nontrivial example is the case 
where T = b/m = const. As a result, we obtain the 
Schrodinger equation with logarithmic nonlinearity [23l | 

ih^ = -^-V 2 x ^+(u + b\n\M 2 y. (IV.37) 

Note that the constant b can be positive, as well as nega- 
tive. From the requirement of the positivity of the inter- 
nal energy, it is easy to obtain the upper bound for the 
case b < 



|6|< 



Anm 2 N 



X. 



(IV.38) 



Concluding this Section, it is worthwhile to reiterate 
that the Schrodinger equation has been derived by the 
sole use of purely classical stochastic arguments. In addi- 
tion, it should be pointed out that the linear Schrodinger 
equation is not the unique possibility. 



V. THE QUANTUM PICTURE 

First of all, we note that g, which has the meaning 
of mass density [see equation <|lV.12[l ] can be rescaled 
g — * mNg, such that it becomes normalized 

/"dxf>(x;t) = l. (V.l) 

This implies that the wave function is normalized as well 

J dx|0(x; t)\ 2 = J dxe(x; t) = 1. (V.2) 



From (|III.23f> and l|IV.15p for the expectation value of an 
arbitrary function -F(x) of position, we obtain 

CF(x)) = ydxdPF(x)/(x,P;t) = ydxF(x)|0(x;i)| 2 . 

The expectation value of momentum can be found in a 
similar manner 

(P) = J dxdPP/(x,P;t) = m J dx(V)g(x;t) 



= -ih J dxi/>*VzV- (V.4) 
Further, we have (for the case, where T = 0) 
(P 2 ) = m 2 J dx£>(x; i)(v 2 + TrC) 



= h- 



J dx(V x 0*) • (V x ip). (V.5) 
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The expectation value for the energy is represented by 
the expression (for the case, where T = 0) 



(H) 



dx 



(V x V) + C/(x;t)|Vr 



(V.6) 



Thus, we recover the basic quantum rules to calculate ex- 
pectation values of observables, which are not higher than 
quadratic in momentum. This implies that not only ex- 
pectation values, but also uncertainties for position and 
momentum can be calculated so that they correspond to 
the standard quantum expression. Similar to the obser- 
vations made by M.J.W. Hall and M. Reginatto [HII3, 
we find that the Heisenberg uncertainty principle can be 
obtained solely by using classical statistical mechanics 
formalism with stochastic ingredients added. 

Since only the mean value of the current velocity v 
and the correlation tensor C mn have been specified to this 
end, it is not immediately clear how one can proceed with 
calculation of expectation values of an arbitrary function 
of momentum. A possible approach to this problem will 
be outlined in the sequel. 

The definition 111.6(1 of the characteristic function can 
be generalized into a characteristic functional 



;t;ix.s:/)) = |^(x;;)| 2 (exp(^s.V 



On the other hand, from 1111.91 we would obtain 



,c7( X .s:/))=V(x+|;t)v*(x-|;t). 



(V.7) 



(V.8) 



analyzed in the Appendix. There, it is shown that equa- 
tion (|V.9|I remains valid to second order in the variable 
s. In addition, an expression for the third-order correla- 
tion function (VkViV,-^ of the random velocity field V is 
derived. 

Let us define the Fourier transform of the wave func- 
tion according to the well-known relation 



^(x;i) 



1 



( , T/ _ r /2_/ dp^(p;t)exp^-x.p). (Y.l.h 



In order to cast a parallel with the discussion in Section 
II, we integrate equation (|V.8(1 over x. As a result, we 
obtain 



G(a)= / dx(0(x,s;t)> 



dp 



cxp [ -s • p 



(Vll) 

According to expressions i|II.8|) and i|II.10|) and their obvi- 
ous generalization, the expectation value of an arbitrary 
function of momentum F(p) can be written in the form 

<F(p)) = F(-ihV s )G(s) = f dpF(p) #(p;f) 2 . 

s=0 J 

(V.12) 

In coordinate representation, we have 



(F(P)) 



dxip*(^;t)F(-iKV x )-ip{ic-,t). (V.13) 



A natural generalization is now in order 



which implies 

|V(x; i)| 2 ^exp ■ V) ^ = ^(x + |; t)f (x - |;i 

(V.9) 

The last equation represents a strong condition to be im- 
posed on the stochastic properties of the random velocity 
field. Clearly, higher order correlation functions of the 
random velocity field V must be specified accordingly in 
order to satisfy equation l|V9(l . The consequences are 



(n*,p)> 



dxV>*(x; t)F(x, -tfiV x )V(x; *), (V.14) 



where as usual, an appropriate operator ordering must 
be specified. 

We end this Section by emphasizing another remark- 
able link between the formalism developed here and the 
Wigner-Weyl approach. Using the compatibility condi- 
tion in the form of (V.8jl . we can rewrite equation 1II.7|I 
as 



J 



(/(x,p;i)) = 



(2nh) n 



dsV>(x+-;*)V>* 



x -H exp H s - p 



(V.15) 



This implies that \ip\ (<5(p — mV)) is the Wigner func- 
tion. Leaving more speculations aside on the fact that 
the latter is a quasi-distribution which is not always pos- 
itive, we note that equation (|V.15J) represents a relation 
between the averaged classical Bohm distribution and the 
Wigner function. 



VI. MOYAL BRACKET 



The Moyal bracket is a useful tool when one wishes to 
determine a semiclassical limit to wave mechanics. Moyal 
|24j elaborated on the theory of Wigner Q on how to de- 
scribe quantum systems in phase space in a way which is 
formally analogous to the dynamics of classical distribu- 
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tions. The Moyal bracket provides a semiclassical limit 
to quantum mechanical commutation relations, which is 
what is of interest to us, and we consider this in some 
detail in the present Section. 

First of all, let us introduce the characteristic dynam- 
ical variable defined as 



C(x, p; k, s) = exp 



-(k • x + s • p) 



(VI.1) 



The rule to calculate expectation values according to ex- 
pression (|V.14|I suggests the introduction of the corre- 
sponding characteristic operator pH ] 



C(x, p; k, s) = exp 



-(k • x + s • p) 



where x implies x, while p = —ih\7 x . The characteristic 
function is defined as the expectation value 



C(k,s) = (C(x,p;k,s)). 



(VI.3) 



According to equation (|V.14(1 and taking into account the 
Campbell-Baker-Hausdorff identity, 



(VI.4) 



when the commutator commutes with both A and B 
(which is the case for operators proportional to x and 
V^), we obtain 



(VI.2) C(k,s) 



dxV>(x+ |;ij<0*(x- | ; t) exp ( ^k • x 



(VI.5) 

Taking into account the inverse Fourier transform of ex- 
pression l|VI.5() and equation (|V.15|) . it immediately fol- 
lows that the characteristic function is a double Fourier 
transform of the Wigner function 



W(x >P ;i) = </(x,p;t)) 



{2irhy 



dkdsC(k, s) exp 



: (k ■ x + s • p) 



(VI.6) 



r 



C(k,s) = / dxdpW(x, p; t) exp 



-(k-x + s-p) 

(VL7) 

Equation i|VI.7|) implies also that the expectation value 
of the characteristic dynamical variable I|V1.1|) is a result 
of integration of its product with the Wigner function 
over all of phase space. 

Further, a generic dynamical variable ^4(x, p) is repre- 
sented by the expression 

i4(x, p) = — / dkds„4(k, s)C(x, p; k, s), (VI. 8) 

(27m) J 



while its corresponding operator in terms of x and p spec- 
ified by the characteristic operator (|VI.2|I can be written 
as 



A(x,p) 



2)1 



dkds^(k,s)C(S,p;k,s). (VL9) 



The dynamical variable l|VI.8|) is usually called a phase 
function of the operator (|VI.9(I . We can take the Fourier 
transform of the phase function ^4(x;p), which returns 
>4(k;s), and substitute it back into equation l|VI.9|) for 
the corresponding operator. This results in 



A(x,p) 



2n 



dxdpdkdsA(x; p) exp 



(x - x) + s • (p - p)] 



(VL10) 



r 



Taking into account expression l|VI.7p for the expectation 
value of the characteristic dynamical variable, we obtain 



(A(x,p)) = (A(%p))= I dxd P A(x,p)^(x,p;t). 

(VI. 11) 



This describes the expectation value of a dynamical 
variable (an operator observable in quantum mechani- 
cal sense) as being the result of integrating its product 
(the product of its corresponding phase function) with 
the Wigner function over all of phase space. In this way 
the Wigner function acts much like a joint probability 
distribution over position and momentum. 
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Expression l|V.14|) shows that the expectation value of which is expected to be proportional to H. We begin by 



the product of two dynamical variables [corresponding 
operators in the sense of equation (|VI.9Jl ]depends on the 
order of the multipliers, and in this sense they do not 
commute. We now wish to determine the measure of non 
commutativity of two arbitrary operators 



ihD(H,p)= A(x,p),B(x,p) , 



(VI. 12) 



J 



substituting in the general expression of equation l|V1.9|) 
for each of the operators A(x, p), B(x,p) and D(x,p) 
This brings us to 



ih / dkdsX>(k,s)C(k,s) = 



2ji 



dkidsidk 2 ds2-4(ki,Si)B(k 2 , s 2 ) C(ki, si),C(k 2) s 2 ) 
I 



(VI. 13) 



where for the sake of simplicity the explicit dependence 
on the operators x and p has been dropped. Taking into 
account the Campbell-Baker-Hausdorff identity (|VI.4|1 . 
and the fact that the commutator 

[ki ■ x + si • p, k 2 ■ x + s 2 • p] = -ih(k 2 ■ si - ki ■ s 2 ), 

(VI. 14) 

is a scalar quantity it is straightforward to verify 



C(ki,si),C(k 2 ,s 2 ) 



We now determine the expectation value on both sides 
of equation (|VI.13(I . resulting in characteristic functions 
on each side of the resulting equation. Substituting ex- 
pression (|VI.7|I for the characteristic function and using 
the representation (|VI.8|I . we obtain 



h / dxdpD(x,p)W(x,p;*) 



2iC(k 1 + k 2 , B X + s 2 ) sin k2 ' Sl m kl ' 82 . (VI.15) 



{2nh) 



^ f dxdpdk 1 ds 1 dk 2 ds 2 T^(x, p; t)A(k u si)B(ka, s 2 ) sin ^ 2 ' Sl ofc kl ' S2 C(k 1; Si)C(k 2 , s 2 ). (VI.16) 



2h 



r 



Next, we note that the argument of the sine on the right- 
hand-side of equation (|VI.16|) can be represented as 

(k 2 ■ si - ki ■ s 2 )C(ki,si)C(k 2 ,s 2 ) 



fi 2 (V xl • V p2 - V p i • Vxa)C(ki,si)C(k2,s 2 ), (VI.17) 



where the subscripts on the differential operators refer to 
which function they operate on, i.e. either C(ki,si) or 
C(k 2 , s 2 ), never both. Bearing in mind that the functions 
«4(ki,si) and S(k 2 ,s 2 ) are independent of x and p, we 
can take the differential operators outside the integrals 
over ki, Si, k 2 and s 2 . As a result, we obtain 



H I dxd P D(x,p)^(x,p;i) = 2 / dxdpW{x,p;t)sm~(V xA • V pS - V pA • V xB )A(x,p)B(x,p), (VI.18) 



where again the subscripts on the differential operators either A(x, p) or B(x, p), never both. Comparing the 
refer to the corresponding function they operate on, i.e. non-Wigner function terms inside the integrals immedi- 
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ately implies Finally, it is worthwhile to mention that the link be- 

tween the formalism used in the deformation quantiza- 
—D(x, p) il0n m °thod and the usual techniques of classical statis- 

2 tical mechanics appears quite natural in our approach. 



= sin-(V Xj4 • V pB - V p a • V x b)A(x, p)B(x, p), 

(VI. 19) 

which is the final result. Considering only the lowest 
order term in h, this reduces to 



£>(x,p) 



dA(x,p) 8B{x,p) dA(x,p) dB{x,p) 



dx 



dp 



dp 



dx ' 
(VI.2 0) 

A common shorthand notation for equation (|VI.19(I is 



APPENDIX A: ANALYSIS OF THE 
COMPATIBILITY CONDITION (TvT9l) 

To verify the validity of the compatibility condition 
(IV.9I . we expand its both sides in a power series in the 
variable s. To third order, we have 



im 

exp ( — s • V 



jD(x,p) = {^(x,p),B(x,p)} Mii , (VI.21) 

where the subscript MB stands for Moyal bracket, as- 
sociating this expression with both the Poisson bracket 
and the commutator bracket. 



VII. DISCUSSION AND CONCLUSIONS 

Starting from a simple classical framework and em- 
ploying some stochastic concepts, the basic ingredients 
characterizing the quantum nature of physical processes 
are recovered. It has been shown that the traditional ax- 
iomatic structure of quantum mechanics can be rebuilt, 
so that the quantum mechanical framework resembles to 
a large extent that of the classical statistical mechanics. 

The main assumption used in the present paper is the 
existence of a random irrotational component in the clas- 
sical momentum. The physical grounds for such assump- 
tion are that an ideal particle localization is not feasible. 
Hence, provided a probability density in configuration 
space is prescribed, the infinitesimal variation of the par- 
ticle "position" in the course of time (i.e., the particle ve- 
locity) cannot be determined precisely. Therefore, there 
is always some uncertainty in the specification of par- 
ticle momentum, which should strongly depend on the 
degree of particle delocalization in configuration space. 
The approach pursued here is by no means an attempt 
to build a realistic model of the underlying momentum 
fluctuations, however some hints concerning their higher- 
order correlation properties are presented. In particular, 
the current velocity fluctuations are shown to be related 
to the turbulent fluctuations in the standard picture of 
Reynolds turbulence. The latter represents an interesting 
and promising guideline for further investigations. 

As a result of the investigation performed, various ba- 
sic elements of the quantum formalism (calculation of 
expectation values, the Heisenberg uncertainty principle, 
the correspondence principle) are recovered by applying 
traditional techniques, borrowed from classical statistical 
mechanics. 



m 2 s k s n n _ im 3 s k sis. 



x (v k viv n + V k Cln + ViC nk + V n C k [ + V k ln) + • ■ ■ , (A.l) 
where 

V kln = (v k ViV n ), (A.2) 

is the yet unknown third order correlation function of the 
random velocity field. Similarly, for the right-hand-side 
of equation l|V.9|) . we obtain 

1, (x + 1 ; t) r (x - 1 ; t) = e + 1 - ) 

+ f^£n^p*Q k Q n1 p _ g k2 p*g n ip _ d n ip*d k ip* + ipd k d n ip*) 



+ fc i n ^ l p*Q k Q l Q n ^ ) _ ^d k did n ip*) 



+ k ' " [dk^did n il)* - d k ip* did n ip + Sym(k,l,n)], 

(A.3) 

where Sym(fc, l,n) denotes four additional terms ob- 
tained by cyclic permutation of the indices fc, / and n. 

Taking into account the explicit expression l|IV.24|) 
for the correlation tensor C kn with T = 0, and the 
amplitude-phase representation 1|III.23|) of the wave func- 
tion, it is a simple matter to verify in a straightforward 
manner that equation (|V.9(I is satisfied up to second or- 
der in the expansion variable s. To satisfy the third order 
however, one needs to specify the third-order correlator 
of the random part of the current velocity. We obtain 
the following expression 

h 2 

v Mn = ? (did„v k + d n d k vi + d k div n ). (A. 4) 
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No attempt is made neither to interpret, nor to provide 
a physical model (which beyond doubt should be nonlo- 
cal) of a possible source, underlying the velocity (momen- 
tum) fluctuations. However, from expression l|A.4|l . it is 
clear that these fluctuations are far from being Gaussian. 
A further insight on the relation between the current ve- 
locity fluctuations and the standard picture of Reynolds 
turbulence is given in Appendix B. 



taking into account equation (|IV.18|) . we obtain 



dC ns 

dt 



dC n 



Vk 



dx k 



dv s 
' dx k 



dv n dV n 



dx k 



1 / dA n OA, 
77i \ at at 



dR 

' dx k 



(B.5) 



APPENDIX B: TURBULENT FLUCTUA TION S 
AND COMPATIBILITY CONDITION (TvT9l) 

To show that the results from Appendix A are consis- 
tent with the model presented in Section IV, we rewrite 
equations (|IV.10|) and (|IV.11|) as 



|f + V x • ( Q V) = 0, 



dV F 

— +V-V,V = -, 

at m 



(B.l) 



(B.2) 



From the assumption that the density g does not fluctu- 
ate, it follows that the fluctuating part of the continuity 
equation HB.ljl reduces to 



(B.3) 



V x • V = -V ■ V X R. 



Averaging the equation for momentum balance (|B.2|) 
and taking into account relation (|B.3|I . we readily obtain 
equation (|1V.18|) . 
Let us further write 

o o -I 

^-(VnVs) + V k — (V n V s ) = -(F n V s + F s V n ). (B.4) 

Ot OXk 777 

which follows directly from the equation for momentum 
balance (|B.2|) . Averaging the last equation l|B.4|) and 



Having already determined the correlation tensor C ns in 
the form ljIV.24Jl . we can manipulate the first term on 
the left-hand-side of equation IB. 5(1 using the continuity 
equation l(IV.17(l . We again consider the case, where T = 
0. The result is 



C s 



—a 



dv n _ dv k 
dx k dx n 



dv s 
dx k 



d 3 



Vk 



d 2 v k dR 

dx n dx s dxk dx n dx s dxk 



dvk 
dx* 



dV n 



dx k 



dR 

dx k 



Vs 



OK 
dt 



V, 



dt 



(B.6) 



If the current velocity v is irrotational as specified by 
(|IV.34|) , the first two terms on the left-hand-side of equa- 
tion (|B.6|) vanish. Assuming also that the current veloc- 
ity fluctuation V is uncorrelated with the random part 
of the force F, we readily arrive (after symmetrization) 
at the expression 1A.4|) for the triple correlation tensor 

'Dnsk • 

Using equations (|B.2I) and l|B.4(l one can proceed in 
calculation of the fourth-order correlator of the current 
velocity fluctuation V. This procedure can be continued 
further and all higher order correlators can be found in 
principle. 
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